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Spectral Structure of Critical Opalescence: Binary Mixture 
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The linearized hydrodynamic equations of irreversible thermodynamics are used to obtain the time 
dependence of the Ath Fourier component of a concentration fluctuation in the one-phase region of 
a binary mixture. The frequency of the light scattered by concentration fluctuations near the critical 
mixing point is found to be broadened into a Lorentzian distribution with a half-width proportional to 
the mass diffusion coefficient, D. The possibility of detecting space dispersion in D is discussed. 

Key Words: Binary diffusion coefficient, concentration fluctuations, critical opalescence, gas 
laser, optical hetrodyne system, spectral distribution of scattered light. 



1. Introduction 

According to the Einstein [1] ^Smoluchowski [2] 
theory of light scattering in a fluid, scattering results 
from fluctuations in the dielectric constant in small 
volumes of the fluid. Since static inhomogeneities 
do not exist in fluids, these fluctuations vary in time 
and this results in a spectrum of frequencies of the 
scattered light [3]. 

Critical opalescence refers to the anomalously large 
amount of scattering which occurs in the vicinity of 
a critical point [4]. In this paper we are concerned 
with the spectrum of critical opalescence in the one 
phase region of a binary mixture. Here the fluctuating 
quantity is the concentration of one of the components 
of the mixture [5]. Large fluctuations in the concen- 
tration are possible near the critical mixing point 
because the chemical potentials of the components 
vary quite slowly with concentration in that region [6], 

Although the spectrum of the scattered light is 
confined to a narrow range of frequencies, it is now 
possible to investigate the spectrum using a gas laser 
and an optical heterodyne system [7]. This type of 
system has been used to examine critical opalescence 
in a cyclohexane and aniline mixture [8]. 

2. Formulation of the Problem 

We shall determine the time dependence of concen- 
tration fluctuations by calculating the decay of a 
fluctuation in time using the linearized equations of 
irreversible thermodynamics [9]. For a binary mix- 
ture, these are the continuity equation 



the Navier-Stokes equation (only the longitudinal 
part of v is included) 



dp/dt + po div v=0, 



(1) 



podv/dt=— grad P + (^ Vs^Vb) grad div v, (2) 

the energy transport equation 

poCpd Tjdt — pokr(d pldc)p yT dc/dt 

+ poT (dSldp) T , cdpld i = kV 2 T, (3) 

and the equation of continuity for the concentration 

podc/d t = poD [ V 2 c + (hrlTo) V 2 T + (k P lp ) V 2 p]. (4) 

In these equations, p is the density, T is the tempera- 
ture, p is the pressure and c is the concentration; 
equilibrium values are denoted by a subscript zero (p ). 
Other quantities are the thermal conductivity k, the 
shear and bulk viscosities r) s and tjb, the diffusion 
coefficient Z), the thermal diffusion ratio kr and the 
barodiffusion ratio k v . The thermodynamic derivatives 
involve the entropy S and the chemical potential p. of 
the mixture. The notation is that of reference 9. 

A light scattering experiment provides a Fourier 
analysis of the fluctuations in the scattering medium. 
If both the scattering angle and the frequency of the 
scattered light are measured precisely, both the space 
and the time variation are Fourier analyzed [3]. Using 
the hydrodynamic equations, eqs (1 to 4), we can deter- 
mine the time dependence of the Ath Fourier compo- 
nent, c(/c, t), of a concentration fluctuation in terms 
of an initial fluctuation c(k). The intensity of the 
scattered light is proportional to (c(k, (o)c( — k)) where 
the (...) indicate an equilibrium average over initial 
(£ = 0) values and 
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/: 



c(k, w) = 2 Re c(k, t)e iw ' dt 



(5) 
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In eq (5) k = 2k x sin 6/2 is the change in the wave vector 
which occurs when light with wave vector fq in the 
medium is scattered through an angle 6. The change 
in the angular frequency is denoted by co. 

It is a straightforward, if involved, procedure using 
eqs (1 to 4) to obtain c(k, t) in terms of the initial fluc- 
tuation c(k) [10]. As expected, two propagating 
(phonon) modes and two diffusive modes are obtained. 
The diffusive modes are a combination of thermal 
diffusion and mass diffusion. Near the critical mixing 
point, where c p is anomalously large and D is small [11], 
the situation is appreciably simpler. There eqs (3) 
and (4) may be used to describe the diffusive modes. 
The pressure terms may be neglected because the 
pressure is insensitive to the local concentration in 
the vicinity of the critical point [6]. 

In the next section, we shall obtain an expression 
for the spectral structure of critical opalescence using 
the constant pressure versions of eqs (3) and (4). 
The results are identical to those obtained using the 
full set of equations. 



3. Calculation 

The equations governing the temperature and the 
concentration are 

dTldt-k T lcp(d fJL/dc) P , T dcldt = (Klp c p ) V 2 T (6) 



Z = (Dk*lc p To)(dnldc)p,T. 
From this it follows that 

c(k, t) _ \9±W-D)+M 



dk) 



-i9±¥x-m±M 



exp[-Wg-<$)k 2 t] 



exp [- \U+<$)kH\ (12) 



where 



6}>=[(x-Dr+2(x+D)t+er- 



-X + D + Z. 



Near the critical point both x an d D(x > D) are small 
[11] so that 



@~ X -D 



and eq (12) simplifies to 



c(k, t) 
c(h) 



= e -Dm + <pm2\ 



(13) 



(14) 



Equation (14) is based on the reasonable assumption 
that (x~ D) 2 vanishes more slowly than D 2 /c p . Com- 
bining eqs (5) and (14) we find that 



dc/dt = D[ V 2 c+(k T /To) V 2 T\ 



(7) 



The first step is to obtain the Fourier (space) and 
Laplace (time) transforms of eqs (6) and (7). The 
transform of the concentration is 



c(k, s) z 



IX 



-ik ■ r„-st 



e-~ st c(r, t)drdt. 



(8) 



Equations (6) and (7) become 

sT(k, s) — sk T /c p (d fJLldc)p, T c(k, s) + (Kk 2 /p c p )T{k, s) 

= T(k) - k r lc p (d fJildc)p, T c(k) (9) 

and 

sc(k, s) + Dk 2 c(k, s) + (Dk T k 2 ITo)T(k, s) = c(k). (10) 

Neglecting T(k) because ( AjTAc) =0 [12], we obtain 



s + x^ 2 + g£ 2 



fhere 



c{k, s) _ 
c(k) " s 2 + skHx + D + i) + xDk 4 



X = K/poCp 



(11) 



2Dk 2 
( dk, co)c( - k) ) = ( c(k)c( - k) ) r (15) 



An equation equivalent to eq (14) has been derived 
by Debye [14J. His derivation is based on the assump- 
tion that there is no dynamical coupling between the 
concentration and the temperature. We have shown 
this to be the case when D < x and both quantities 
are small. Then it is possible to separate the two 
diffusive modes into one dominated by the thermal 
diffusity x an d one dominated by the diffusion coeffi- 
cient D. The mode dominated by the thermal dif- 
fusity has effectively zero weight in the sum rule, eq 
(16), for the intensity of the scattered light. 



\ - f °° < c(&, co)c( - k) )da> = ( c(k)c( - k) ) . 

2 TT J- x 



(16) 



(The phonon modes, which we have not explicitly 
considered, also have zero weight.) 

If only the intensity of the scattered light is observed, 
( c(k)c( — k) ) is obtained. Near the critical point 
( c(k)c( — k) ) is commonly assumed to be of the form 
(A- 2 + k 2 ) -1 where k 1 is the range of the two particle 
correlation function [13]. Away from the critical point, 
( c(k)c( — k) ) is independent of the wave vector and 
equals the mean square fluctuation in the concen- 
tration. 
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4. Discussion 

The spectrum of critical opalescence is that of a 
Lorentzian distribution with a half width of Dk 2 /n 
cycles/sec. Recent measurements of the spectrum 
of critical opalescence in a mixture of cyclohexane 
and aniline are of this type [8|. Measurements of the 
spectral structure of critical opalescence provide 
a gradient-free method of determining the diffusion 
coefficient. 

The type of calculation presented above is certainly 
appropriate when the range of the two particle cor- 
relations is small compared to the wavelength of the 
light which is scattered. Near the critical point, the 
range of the correlations is on the order of the wave- 
length of visible light. It is not unreasonable to expect 
that near the critical point, transport coefficients 
may reflect the long range of these correlations by 
becoming /-dependent [15]. This interesting pos- 
sibility could be checked by seeing whether or not 
the half width of the Lorentzian distribution varies 
linearly with fr. Any deviation from linearity would 
suggest that the diffusion coefficient I) should be re- 
placed by D{k). a /.-dependent diffusion coefficient. 
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